Direct Measurement of the Fermi Energy in Graphene Using a Double Layer 

Structure 
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We describe a technique which aUows a direct measurement of the relative Fermi energy in an 
electron system using a double layer structure, where graphene is one of the two layers. We illustrate 
this method by probing the Fermi energy as a function of density in a graphene monolayer, at zero 
and in high magnetic fields. This technique allows us to determine the Fermi velocity. Landau level 
spacing, and Landau level broadening in graphene. We find that the A'^ = Landau level broadening 
is larger by comparison to the broadening of upper and lower Landau levels. 

PACS numbers: 73.43.-f, 71.35.-y, 73.22.Gk 



The Fermi energy is a fundamental property of an elec- 
tron system, and thermodynamic measurements which 
probe the Fermi energy or density of states are key to 
understanding the host material band structure and elec- 
tron interaction effects. Although a number of thermo- 
dynamic properties, such as specific heat [Ull], magneti- 
zation [3], magnetocapacitance or compressibility [11 
can directly probe the density of states in an electron 
system, accessing them experimentally becomes increas- 
ingly difficult at the micro- and nano-scale. In the case 
of graphene 0] , magnetization and specific heat measure- 
ments are exceedingly difhcult, and the accuracy of com- 
pressibility [vj and capacitance measurements |8l-[lo| are 
also limited by the reduced sample dimensions. Using a 
double layer device structure where graphene is one of 
the layers, we describe a technique which allows a direct 
measurement of the Fermi energy in an electron system 
with an accuracy which is independent of the sample size. 
The underlying principle of the method discussed here is 
that an interlayer bias applied to bring the graphene layer 
to the charge neutrality point is equal to the Fermi en- 
ergy of the opposite electron system. We illustrate this 
technique by probing the Fermi energy in a graphene 
layer, both at zero and in high magnetic fields. We show 
that this method allows an accurate determination of the 
Fermi velocity in graphene, the Landau level spacing, and 
Landau level broadening in high magnetic fields. 

Our samples are independently contacted graphene 
double layers, consisting of two graphene single layers 
separated by a thin dielectric as shown in Fig. 1(a) 11 1. 
To fabricate such devices, we first mechanically exfoliate 
the bottom graphene layer from natural graphite onto a 
280 nm thick Si02 dielectric, thermally grown on a highly 
doped Si substrate. Standard e-beam lithography, Cr/Au 
deposition followed by lift-off, and O2 plasma etching are 
used to define a Hall bar device. A 4 to 7 nm top AI2O3 
dielectric layer is deposited on the bottom layer by atomic 
layer deposition, and using evaporated Al as a nucleation 
layer. The dielectric film thickness grown on graphene is 
further verified by transmission electron microscopy in 



multiple samples. To fabricate the graphene top layer, 
a separate graphene single layer is mechanically exfoli- 
ated on a Si02/Si substrate. After spin-coating poly- 
metyl metacrylate (PMMA) on the top layer and curing, 
we etch the underlying substrate with NaOH, and detach 
the top layer along with the alignment markers captured 
in the PMMA membrane. The membrane is transferred 
onto the bottom layer device, and aligned. A Hall bar 
is subsequently defined on the top layer, completing the 
double layer graphene device. Three samples were inves- 
tigated in this study, all with similar results. We focus 
here on data collected from one sample with a 7.5 nm 
thick interlayer dielectric, and with an interlayer resis- 
tance larger than 1 Both layer mobilities are 10,000 
cm^/V-s. Using small signal, low frequency lock-in tech- 
niques we probe the layer resistivities as a function of 
back-gate bias {Vbg)i and the inter-layer {Vtl) bias ap- 
plied on the top layer. The bottom layer is maintained 
grounded during measurements. 

Figure l(b,c) data show the longitudinal resistivity of 
the bottom {pb) and top {pr) layer measured as a func- 



tion of Vtl, and at different Vbg values [12|. The data 
Pb,t vs. Vtl exhibit the ambipolar behavior charac- 
teristic of graphene, and with a charge neutrality point 
which is VsG-dependent. The shift of the charge neu- 
trality point of the bottom layer as a function of Vbg is 
explained by picturing the bottom layer as a dual-gated 
graphene single layer, with the Si substrate as back-gate 
and the top graphene layer serving as top-gate. The de- 
pendence of the pt vs. Vtl data on Vbg is more subtle, 
and implies an incomplete screening by the bottom layer 
of the back-gate induced electric field. 

We can quantitatively explain the top and bottom 
layer density dependence on Vbg and Vtl using a simple 
band diagram model. Figure 1(d) shows the band dia- 
gram of the graphene bilayer at a finite VbGi and with 
both layers at ground potential. For simplicity the back- 
gate Fermi energy and the two graphene layers charge 
neutrality points are assumed to be at the same energy at 
Vbg = V. Once a finite Vbg is apphed, charge densities 



2 



(a) v.. 




- AI203 

- SiOj 
-Si 



irJfmJi 




Back I Bottom Top ©eVs/O; 
Gate Graphiene Graphere (Df^.^ 



FIG. 1: (color online) (a) Schematic representation of a 
graphene double layer separated by an AI2O3 dielectric, and 
with a bottom Si02 dielectric. A back-gate (Vbg) and inter- 
layer (Vtl) bias can be applied on the Si substrate and top 
layer, respectively. Lower right: optical micrograph of a 
graphene double layer device. The red (blue) contour marks 
the bottom (top) layer. The scale bar is 5 /xm. (b,c) Layer re- 
sistivities measured as a function of Vtl and Vbg at T = 0.4 
K. (d,e) Band diagram of a graphene double layer under an 
applied back-gate [panel (d)] or inter-layer [panel (e)] bias. 



are induced in both bottom {ns), and top (jt-t) layers. 
Consequently electric fields are built-in across both bot- 
tom Si02 and inter-layer AI2O3 dielectric. The applied 
Vbg potential is the sum of the potential drop across the 
Si02 dielectric and the Fermi energy of the bottom layer: 



eVsG = e^{nB + nT)/Csi02 + EpiriB) 



(1) 



Here, Ep{nB) represents the Fermi energy of the bot- 
tom layer corresponding to a charge density n b , and mea- 
sured from the charge neutrality point; nB,T and Ep(n) 
are positive when the carriers are electrons, and negative 
when the carriers are holes. CsiO^ denotes the bottom 
dielectric capacitance per unit area. 



Figure 1(e) shows a similar band diagram of the 
graphene bilayer, but in the presence of a finite inter- 
layer bias and at Vbg = V. Similarly to Fig. 1(d), 
the applied Vtl bias can be written as the sum of the 
potential drop across the AI2O3 dielectric and the Fermi 
energies of the two layers: 



cVtl 



= EpiuB) - [EFinr) + e^nT/CAi^o,) (2) 

Here, Ep{tit) represents the Fermi energy of the top 
layer at a charge density ut, and CAI2O3 is the interlayer 
dielectric capacitance per unit area. In Fig. 1(e) the 
Vtl bias is assumed to be positive, resulting in electrons 
(holes) induced in the bottom (top) layer. Although we 
derived Eqs. (1) and (2) assuming Vtl = V, and 
Vbg = V respectively, the two equations hold at all 
Vbg and Vtl values. Most importantly, we do not make 
any assumption with regard to the Ep dependence on 
ub and jit- As we show below, this dependence will be 
determined experimentally. 

Figure 2(a) data show contour graphs of pB (left panel) 
and pt (right panel) as a function of Vbg and Vtl- The 
bottom layer resistivity dependence on gate bias is very 
similar to a dual-gated graphene monolayer ll|, show- 



ing an almost linear dependence of the charge neutral- 
ity point on Vbg and Vtl, with a slope equal to the 
CsiO^/CAhOs ratio. Using Csio^ = 12 nF/cm^ for the 
bottom Si02 dielectric, we determine the inter-layer di- 
electric capacitance to be CAhO-i = 340 nF/cm'^. The 
capacitance values are confirmed by Hall measurements. 

The top layer resistivity shows the characteristic am- 
bipolar behavior as a function of Vtl, and with a weaker 
Vbg dependence. Let us examine more closely the top 
layer charge neutrality point dependence on Vbg and 
Vtl- If we consider the top layer at the charge neutrality 
point, setting riT = in Eq. (2) yields: 



cVtl = EF{nB) 



(3) 



This equation contains a simple, yet remarkable result. 
The inter-layer bias required to bring the top layer at the 
charge neutrality point is equal to the Fermi energy of the 
opposite layer, in units of eV [Fig. 2(b)]. Consequently, 
tracking the top layer charge neutrality point in the Vbg 
- Vtl plane [dash-dotted trace in Fig. 2(a) left panel], 
results in a measurement of the bottom layer Fermi en- 
ergy as a function of Vbg- Furthermore, setting tlt — 
in Eq. (1), and using Eq. (3) allows for to be deter- 
mined as a function of Vbg and Vtl along the top layer 
charge neutrality line of Fig. 2(a): 



Vbg - Vtl = euB/CsiO^i 



(4) 



Equations (3) and (4) provide a direct measurement of 
the bottom layer Fermi energy as a function of density. 
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FIG. 2: (color online) (a) Contour plots of ps (left panel) 
and pT (right panel) measured as a function of Vbg and Vtl, 
at T = 0.4 K. The bottom layer responds to Vbg and Vtl 
similar to a dual-gated monolayer, (b) Band diagram of a 
graphene double layer with the top layer at the charge neu- 
trality point. The inter-layer bias is equal to the bottom layer 
Fermi energy, (c) Bottom layer Fermi energy vs. carrier con- 
centration, determined using data of panel (a). The symbols 
are experimental data, while the solid line represents the ex- 
pected fiVFy/TinB dependence. 



To illustrate this, in Fig. 2(c) we show the bottom layer 
Fermi energy Ep B as a function of n^, determined using 
Fig. 2(a) data and Eqs. (3) and (4). The Ep values are 
in excellent agreement with the Ep^ns) — Tivp^imB de- 
pendence expected for the linear energy-momentum dis- 
persion of graphene, and with an extracted Fermi velocity 
of vp — 1.15 X 10^ cm/s. 

In the following we show that the above method applies 
equally well to an electron system in high magnetic fields, 
allowing a direct measurement of Landau level (LL) en- 
ergies and broadening. In Fig. 3(a) we show the contour 
plots of pt (top panel) and pB (bottom panel) measured 
as a function of Vbg and Vtl in an applied perpendic- 
ular magnetic field i? = 8 T. Both layers show quantum 
Hall states (QHS) marked by vanishing resistitvities at 
filling factors v = A{N + \), consistent with a graphene 
monolayer [l^, [13] • The integer N represents the Lan- 
dau level index. The top panel of Fig. 3(a) data shows 
a step-like dependence of the top layer charge neutrality 
point on Vbg and Vtl- Similarly to Fig. 2, substituting 
cVtl with Ep^B at the top layer charge neutrality line in 
Fig. 3(a) (top panel) provides a mapping of Ef^b as a 
function of Vbg- To visualize this, the top layer charge 
neutrality line in the Vbg — Vtl plane is superposed with 
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FIG. 3: (color online) (a) pr (top) and pB (bottom) con- 
tour plots measured as a function of Vbg and Vtl at B = 8 
T, and T = 0.4 K. Both layers show quantum Hall states 
marked by vanishing longitudinal resistance at filling factors 
V — ±2, 6, 10, consistent with mono-layer graphene. The top 
layer charge neutrality line (dashed line) shows a step-like be- 
havior, with the steps matching the bottom layer quantum 
Hall states, (b) pB (blue line, top axis) vs. Ef,b = sVtl, 
and Ef,b vs. ub (red line, bottom axis) determined from the 
top layer charge neutrality line of panel (a). The Ef,b values 
at the peak positions of ps provide the Landau level energies. 



the Pb contour plot of Fig. 3(b) (bottom panel), which 
shows step-like increments of Ep^B coinciding with the 
QHSs of the bottom layer. 

Figure 3(b) shows Ep,B vs. ris at -B = 8 T deter- 
mined by tracking the top layer charge neutrality line in 
the Vtl — Vbg plane in Fig. 3(a), and using Eqs. (3) and 
(4) to convert Vtl and Vbg into Ep^B and ub, respec- 
tively. In addition. Fig. 3(b) shows pB vs. Ep^B, deter- 
mined by tracking the bottom layer resistivity along the 
top layer charge neutrality line [dashed-dotted line of Fig. 
3(a)]. Figure 3(b) data manifestly shows the staircase- 
like behavior expected for the Fermi level dependence on 
density for a two-dimensional electron system in a per- 
pendicular magnetic field. The peaks in the pB vs. Ep^B 
data of Fig. 3(b), corresponding to the Fermi level lying 
in the LL center and probing extended states, correlate 
with plateaus in the Ep ^ vs. ub, associated with the 
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large LL density of states. The peaks in the ps vs. Ef,b 
data of Fig. 3(b) provide a direct measurement of the 
LL energy. Figure 4(a) summarizes the bottom graphene 
layer LL energy as a function of index {N) at B = 8 
T. The experimental data is in excellent agreement with 
the theoretical dependence — ±.VF\/'2.fieB\N\, cor- 
responding to a Fermi velocity vp = 1.17 x 10® cm/s, 
a value less than 2% different than the Fermi velocity 
determined at i? = T using Fig. 2 data. 

In Figure 4(b) we compare the Ef.b vs. ub data deter- 
mined experimentally at i? = 8 T, with calculations. As- 
suming a Lorentzian distribution of the disorder-induced 
LL broadening, the density of states D{E) writes: 



D(E) 



h ^ ' TT 



N 



TT {E- 



(5) 



with being the broadening of the A^-th LL. The carrier 
density (n) dependence on Ep in the limit T = K is: 



niEp) = 



D{E)dE 



(6) 



Using Eqs. (5) and (6), the best fit to Ef.b vs. nB 
data is obtained for 70 = 14 meV, and — 6.5 meV 
for |7V| > 0. The summation in Eq. (5) does not con- 
verge if carried out to infinity, and a high-energy cut-off 
is customarily used. For the calculations of Fig. 4(b) 
we used \N\ < 100 in Eq. (5), corresponding to a 1 eV 
cut-off energy; increasing the cut-off LL index to 1,000, 
will change the best fit 7 value by less than 0.5 meV. 
The lower inset of Fig. 4(b) shows a comparison of the 
Ef,b vs. ng experimental data with calculations using 
the same broadening for the = LL as the upper and 
lower LLs, 7 = 6.5 meV. The larger broadening of the 
A^ = LL by comparison to the other LLs is an interest- 
ing finding. A theoretical study 15| , which examined the 
impact of static disorder on LL broadening in graphene 
without considering interaction showed that the A^ = 
LL broadening is the same as for the other LLs. On 
the other hand electron-electron interaction can impact 
the broadening of the four-fold degenerate A^ = LL, 
and experimental data on exfoliated graphene on Si02 
substrates show a splitting of the A^ = LL in high, 
B = 45 T magnetic fields 16| , explained as a many-body 
effect. Lastly, we note that a Gaussian-shaped Landau 
levels density of states yields worse fits to Fig. 4 data, 
by comparison to the Lorentzian shape density of states. 
Scanning tunneling microscopy studies [13, [Si , and com- 
pressibility measurements in graphene [7| also favor the 
Lorentzian LL lineshape by comparison to the Gaussian 
one. A recent theoretical study argues that LL local den- 
sity of states has a Lorentzian lineshape while the total 
density of states is Gaussian [l^. Presumably, the sam- 
ple size examined here, defined by a 4 /im Hall bar width 
coupled with the 8 ^m top layer contact spacing is suf- 
ficiently small such that the Lorentzian LL line-shape 
dominates. 
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FIG. 4: (color online) (a) Landau level energy in mono-layer 
as a function of index (N). The symbols are experimental 
data determined from the Ef,b positions at the ps peaks in 
Fig. 3(b). The solid line is the theoretical ±VFy^2heB\N\ 
dependence corresponding to vf = 1.17 x 10* cm/s. (b) Bot- 
tom layer Fermi energy vs. density at _B = 8 T. The symbols 
represent experimental data, and the solid (red) line is a fit 
assuming a Landau level Lorentzian line shape. The best fit 
is obtained for a LL broadening of 7jv = 6.5 meV for | A| > 0, 
and 70 = 14 meV. The upper inset shows the calculated den- 
sity of states corresponding to the best fit to experimental 
data. The lower inset shows the Ep vs. ub data in the vicin- 
ity of zero density. The symbols represent experimental data, 
and the lines are calculations assuming 70 — 6.5 meV (blue 
line) and 70 — 14 meV (red line). 



In summary, we present a method to determine the 
Fermi energy in a two-dimensional electron system, using 
a double layer device structure with one layer consisting 
of graphene. We illustrate this technique by probing the 
Fermi energy in a separate graphene layer as a function 
of density at zero, and in a high magnetic field, and de- 
termine with high accuracy the Fermi velocity, and the 
Landau level broadening. The technique sensitivity is 
independent of the sample dimensions, which makes it 
applicable to a variety of nanoscale materials. 
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